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Creating Two-Dimensional Electron Gas in Nonpolar Oxide Interface via Polarization 
Discontinuity: First-Principles Analysis of CaZrOs/SrTiOs Heterostructure 
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We studied strain-induced polarization and resulting conductivity in the nonpolar/nonpolar 
CaZrOs/SrTiOs (CZO/STO) heterostructure (HS) system by means of first-principles electronic 
structure calculations. By modeling four types of CZO/STO HS-based slab systems, i.e., Ti02/Ca0 
and Sr0/Zr02 interface models with CaO and Zr02 surface terminations in each model separately, 
we found that the lattice-mismatch-induced compressive strain leads to a strong polarization in the 
CZO film, and as the CZO film thickness increases, there exist an insulator-to-metal transition. The 
polarization direction and critical thickness of the CZO film for forming interfacial metallic states 
depend on the surface termination of CZO film in both types of interface models. In the Ti02 /CaO 
and Sr0/Zr02 interface models with CaO surface termination, the strong polarization drives the 
charge transfer from the CZO film to the first few Ti02 layers in the STO substrate, leading to the 
formation of Two-Dimensional Electron Gas (2DEG) at the interface. In the HS models with Zr02 
surface termination, two polarization domains with opposite directions were formed in the CZO 
film, which results in the charge transfer from the middle CZO layer to the interface and surface, 
respectively, leading to the co-existence of the 2DEG on the interface and the Two-Dimensional 
Hole Gas (2DHG) at the middle CZO layer. These findings open a new avenue to achieve 2DEG 
(2DHG) in perovskite-based HS systems via polarization discontinuity. 

PACS numbers: 73.20.-r, 73.40.-c, 71.15.-m 


I. INTRODUCTION 

The perovskite-based oxide heterostructures (HS) are 
attracting increasing interests because of their novel in¬ 
terfacial properties such as the interfacial superconduc¬ 
tivity and ferromagnetism that are drastically differ¬ 
ent from those of the corresponding bulk materials. [1- 
4] One typical example is the formation of the 
high-mobility Two-Dimensional Electron Gases (2DEG) 
at Ti02-terminated interface in the polar/nonpolar 
LaA 103 /SrTi 03 (LA 0 /ST 0 ) HS system. [5, 6] The ex¬ 
cellent conducting property of the 2DEG at the inter¬ 
face of the HS makes it promising for applications in the 
next-generation nanoelectronic oxide devices. [7—10] Gen¬ 
erally speaking, there are three main approaches to form 
a 2DEG in a HS system. The first one is via modu¬ 
lation doping, [11] which often creates a 2DEG in a HS 
that consists of a wide-band-gap and a narrow-band-gap 
semiconductor. The n-type doping in the wide-band-gap 
semiconductors introduces additional electrons in the HS, 
and the electrons transfer from the high conduction band 
edge of the wide-band-gap semiconductor to the conduc¬ 
tion band of the narrow-band-gap semiconductor on the 
other side of the HS leads to the 2DEG.[iO] One exam¬ 
ple of the 2DEG utilizing this mechanism is the band 
bending AlGaAs/GaAs HS system, where the electrons 
forming the 2DEG in the interfacial GaAs layer are sup¬ 
plied from the doped AlGaAs.[12, 13] The second one is 
by means of electronic reconstruction caused by the polar 
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discontinuity at the interface between a polar and a non¬ 
polar layer in a HS system such as the LAO/STO. In the 
LAO/STO HS system, the charge transferred from the 
polar (LaO)^'*' layers to non-polar (Ti02)° layers par¬ 
tially occupy the Ti 3d orbitals in the (Ti02)° layers 
near the interfacial region, thus forming the 2DEG in 
the STO layer. [14-16] The third one is via a disconti¬ 
nuity of the spontaneous or strain-induced polarization 
that creates an internal electrostatic field and further pro¬ 
duces polarization charges bound at the interface, creat¬ 
ing the 2DEG. One typical example 2DEG system using 
this mechanism is the traditional semiconductor material 
Zni_a;Mg2;0/Zn0. [17-19] In the Zni_a;Mga;0/ZnO HS 
system, the polarization-induced internal electrostatic 
filed produces and confines the uncompensated bound 
charge at the interface, forming the 2DEG.[20] In addi¬ 
tion to these three mechanisms, several other factors such 
as oxygen vacancies [2 1-25] and cation intermixing[14, 26] 
were also proposed to be able to produce metallic states. 
For instance, it was recently revealed that tuning redox 
reactions on the surface of STO substrate might be one 
approach to produce oxygen vacancies and the resulting 
interfacial conductivity. [21] 

Compared to the success of generating 2DEG in the 
LAO/STO system via the polar discontinuity, [5, 6] there 
have been few reports on the possibility to produce 
the 2DEG in the perovskite oxide HS using the polar¬ 
ization. Very recently, Chen and co-workers prepared 
CaZrOs/SrTiOs (CZO/STO) perovskite oxide HS by de¬ 
positing the CZO film on the (001) Ti02-terminated STO 
substrate using Pulsed Laser Deposition (PLD) and ob¬ 
served an insulator-to-metal transition in the HS system, 
when the CZO film thickness t increases up to around 
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6 unit cells (uc).[27] The CZO/STO HS shows metallic 
behaviour at t > 6 uc but insulating at t < 6 uc. By com¬ 
binatorial Scanning Transmission Electron Microscopy 
(STEM) and Electron Energy-Loss Spectroscopy (EELS) 
characterization, they found that the Ca^’*' and Zr^+ 
cations in the CZO film move towards the STO sub¬ 
strate, which is caused by the lattice-mismatch-induced 
compressive strain on the CZO film. This indicates that 
the relative displacement between the cations and an¬ 
ions in the CZO film produces a piezoelectric polariza¬ 
tion, which plays a crucial role in forming the 2DEG in 
the CZO/STO HS system. This experimental finding, 
for the first time, has provided a solid proof to produce 
a high-mobility (^60,000 cm^ s“^ at 2K) 2DEG in 
the non-polar perovskite oxide HS using the piezoelectric 
polarization effects. In spite of the encouraging experi¬ 
mental finding, several open fundamental questions arise 
naturally such as the influence of the surface termination 
on the interfacial electronic properties and the depen¬ 
dence of the electron transport property on the thickness 
of the CZO film. Therefore, from the theoretical view¬ 
point, it is still very necessary to deeply understand the 
underlying working mechanism to produce the interfacial 
2DEG in the CZO/STO HS system. 

In this work, to reveal the underlying mechanism to 
produce the interfacial conductivity in the CZO/STO HS 
system, we studied the structural and electronic prop¬ 
erties of the HS system using first-principles electronic 
structure calculations. Four types of CZO/STO HS- 
based slab models, i.e., Ti 02 /Ca 0 and Sr0/Zr02 in¬ 
terface models with CaO and Zr02 surface terminations 
in each model, were studied. We investigated the strain- 
induced polarization and resulting electronic properties 
in these models, and found that there exists an insulator- 
to-metal transition as the CZO film thickness increases. 
The relative displacement between the cations and an¬ 
ions, polarization direction and strength, charge transfer, 
and the critical thickness of the CZO film for forming the 
interfacial metallic states were elucidated. This work pro¬ 
vides a clear picture for the strain-induced polarization 
and consequent formation of the 2DEG in the nonpolar 
perovskite-based HS system. 


II. STRUCTURAL MODELING AND 
COMPUTATIONAL DETAILS 

In principle, two types of interface structures, i.e., 
Ti02/Ca0 and Sr0/Zr02, can be modeled. In each 
model, the CZO film is deposited on the Ti02-terminated 
and SrO-terminated (STO)io substrate along the [001] 
direction, respectively, in which 10 uc of the STO was 
used to model the substrate. Considering that differ¬ 
ent surface termination of the CZO film, i.e., the CaO 
surface and the Zr 02 surface, may influence the polar¬ 
ization and the resulting electronic property of the HS 
system, we added a vacuum layer of about 15 A along 
the [001] direction of the HS model above the CZO film. 


Hence, to perform a comprehensive comparison study, 
four types of HS-based slab models, i.e., Ti02/CaO in¬ 
terface models with CaO and Zr02 surface termination, 
and SrO/Zr02 interface models with CaO and Zr02 
surface termination, were simulated in this work. An¬ 
other reason for doing this is based on the experimental 
consideration. Normally, that is, under stoichiometric 
condition, the Ti 02 /CaO interface model is terminated 
with Zr02 surface, and the SrO/Zr02 interface model 
is terminated with CaO surface. Actually, by tuning 
the chemical potential of the composition of the CZO 
film in the experiment, that is, under non-stoichiometric 
condition, the Ti 02 /CaO interface model may be ter¬ 
minated with CaO surface, and similarly, the SrO/Zr02 
interface model may be terminated with Zr 02 surface. 
A recent experiment revealed that, in the well-known 
(Ti02)°/(LaO)^''' LAO/STO HS system, the Sr atoms 
can even diffuse from the interface to the surface, form¬ 
ing SrO surface layer. [28] 

As discussed later, after structural relaxation, strong 
polarization may occur in the CZO/STO HS models. The 
periodic boundary condition in the supercell approach 
determines that there might exist an electric-field inter¬ 
action between the surface of the STO substrate and the 
surface of the CZO film, though the strength strongly de¬ 
pends on the vacuum thickness. To exclude the influence 
of such type of electric-filed interaction on the polariza¬ 
tion and electronic property, we made test calculations by 
inserting a 40 A vacuum layer into the Ti 02 -terminated 
CZO/STO HS model, in which the electrostatic interac¬ 
tion between the neighboring surfaces can be negligible. 
Our test calculations show that the two CZO/STO HS 
models with a 15 A vacuum layer and a 40 A vacuum 
layer give exactly same results, suggesting that a vac¬ 
uum of 15 A in the CZO/STO HS-based slab system is 
appropriate to avoid the presumably electrostatic inter¬ 
action. 

The STO crystal has a cubic phase (space group no. 
221, PmSm) at room temperature with a lattice parame¬ 
ter of 3.905.A[.5] The CZO crystallizes in an orthorhombic 
(space group no. 62, Pnma) at room temperature, [29] 
and there exists a tilting and rotation of the ZrOe oc- 
tahedra (corresponding to Zr-O-Zr bending). [30] To de¬ 
scribe the epitaxial perovskite-based heterostructure, one 
convenient approach is to treat the distorted perovskite 
as pseudo-cubic phase. In this work, the distorted CZO 
crystal was modeled in a pseudo-cubic phase with a lat¬ 
tice constant of 4.012 A, [27], which lead to a lattice mis¬ 
match of approximately 2.7%. To resemble the experi¬ 
mental epitaxial material growth process, all the ions are 
fully relaxed by minimizing the atomic forces up to 0.03 
eV /A, while the cell parameters along the x- and y-axis 
were fixed. 

Our Density Functional Theory (DFT) electronic 
structure calculations were carried out using the Vi¬ 
enna Ab-initio Simulation Package (VASP).[31, 32] The 
Projector Augmented Wave (PAW) potentials were 
used for electron-ion interactions [33] and the General- 
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FIG. 1. (Color online) Calculated total DOS for unrelaxed 
(first column) and relaxed (second column) (CZO)n/STO (n 
= 3.5, 4.5, 5.5, and 6.5) HS-based slab (Ti02/Ca0 interface 
model) with CaO surface termination. The vertical dashed 
line indicates the Fermi level in this and each subsequent DOS 
plot. 


ized Gradient Approximation (GGA) parameterized by 
Perdew-Burke-Ernzerhof (PBE) was used for exchange- 
correlation functional. [34] The cut-off energy of 450 eV 
for the plane wave basis set and the 10 x 10 x 1 fc-space 
grid in the irreducible wedge of the Brillouin zone were 
used. The self-consistency was assumed for a total en¬ 
ergy convergence of less than 10“® eV. The density of 
states (DOS) was calculated using the Gaussian smear¬ 
ing of 0.05 eV. 


III. RESULTS AND DISCUSSION 
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FIG. 2. (Color online) Calculated (a) layer-resolved partial 
DOS for relaxed (CZO) 6 . 5 /STO slab (Ti02/CaO interface 
model) with CaO surface termination, along with the charge 
density projected on the bands forming the metallic states 
near the interface and surface region, (b) Orbital-resolved 
partial DOS for Ti atoms at the IF-I Ti02 layer and for O 
atom at the surface CaO layer. 


A. Ti02/Ca0 Interface Model 

1. CaO Surface Termination 

We first examined the polarization and electronic prop¬ 
erty of the Ti02/CaO interface model with CaO surface 
termination. To investigate the critical thickness of the 
CZO him for forming interfacial metallic states, we mod¬ 
eled the CZO/STO HS system by depositing the CZO 
him with various uc on the Ti02-terminated STO sub¬ 
strate along the [001] direction. Hereafter, these interface 
systems are referred to as (CZO)„/STO, in which n de¬ 
notes the number of the CZO uc. To explicitly show the 
polarization-induced electronic property modihcation, we 
calculated the total DOS for unrelaxed (hrst column) 
and relaxed (second column) (CZO)„/STO (n = 3.5, 4.5, 
5.5, and 6.5) HS-based slab system in Fig. 1. Our cal¬ 


culated total DOS clearly exhibits an insulating nature 
for all the unrelaxed CZO/STO HS-based slab systems. 
This implies that there is no charge transfer between the 
nonpolar CZO him and STO substrate, unlike the po¬ 
lar/nonpolar LAO/STO HS system. In contrast, the re¬ 
laxed CZO/STO HS models (second column in Fig. 1) 
display an insulator-to-metal transition when the num¬ 
ber of the CZO uc n increases up to 6.5. For the models 
with the CZO uc n less than 6.5, all the systems are in¬ 
sulating. This indicates that the critical thickness of the 
CZO him for forming the conducting state is about 6.5 
CZO uc, which is in an excellent agreement with the re¬ 
cent experimental hnding that the CZO/STO HS model 
shows metallic behavior at u > 6 uc.[27] 

To reveal the origin of metallic states in the CZO/STO 
HS system, we calculated layer-resolved DOS of the 
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(CZ0)6.5/ST0 HS system with CaO surface termina¬ 
tion in Fig. 2a. For a direct view of the each layer’s 
contribution to the metallic states, we also plotted the 
three-dimensional (3D) charge density projected on the 
bands forming the metallic states near the Fermi level. 
Our results show that O 2p states in the CZO layers sig¬ 
nificantly shift toward higher energy when the CZO layer 
moves from the interface to the vacuum. At the CaO sur¬ 
face layer, occupied O 2p states cross the Fermi level, pro¬ 
ducing some unoccupied states. The O 2p states from the 
2nd Zr02 and 3rd CaO surface layers also cross the Fermi 
level by producing minor unoccupied states. This indi¬ 
cates that the CZO/STO HS forms p-type conducting 
states on the surface. In contrast, at the interface, some 
Ti 3d states cross the Fermi level by forming some oc¬ 
cupied states, thus producing n-type conducting states. 
For convenience of the discussion, the 1st, 3rd, and 5th 
layers of STO was defined as IF-I, IF-III, and IF-V, re¬ 
spectively. Our calculations clearly show that most of the 
interfacial metallic states near the Fermi level come from 
the IF-I Ti02 layer of the STO substrate, along with a 
small contribution from the IF-III and IF-V Ti02 lay¬ 
ers. Further partial DOS analysis indicates that all the 
other layers away from the interface exhibit an insulating 
behavior, showing a typical character of the 2DEG. 

The calculated charge density plot in Fig. 2a clearly re¬ 
flects the same pattern that the conducting states mostly 
come from three sequential Ti 02 layers close to the inter¬ 
face, surface CaO layer, and subsurface Zr02 and CaO 
layers. The geometrical structure in the charge density 
plot shows that the compressive strain exerted on the 
CZO film leads to the downward movement of the Ca^+ 
and Zr^+ cations towards the interface, causing a rela¬ 
tive displacement between the cations and anions, thus 
producing a strong “down” polarization, Pdn- A slight 
downward distortion of the Ti cation away from the in¬ 
terface also occurs in the first few layers of the STO 
substrate, and thus there exists a polarization discon¬ 
tinuity between the CZO film and STO substrate. As a 
result, the strong polarization produces an internal elec¬ 
trical field in the CZO film, which leads to the charge 
transfer from the occupied O 2p orbitals in the surface 
to the Ti 3c? orbitals at the interface, though the whole 
CZO/STO HS system keeps the charge neutrality. The 
transferred electrons are mainly captured by the local¬ 
ized Ti 3d orbitals of the first three Ti02 layer in the 
STO substrate, which forms n-type interface conducting 
states. Moreover, our calculations show that the Zr 4c? 
orbital substantially shifts towards the lower energy as 
the the Zr 02 layers move from the surface to the inter¬ 
face (a line is marked in the STO layers for a guide to 
the eye), though it has no contribution to the interfa¬ 
cial conductivity. This implies a significant band bend¬ 
ing between the STO substrate and the CZO film. In 
addition, the calculated orbital-resolved partial DOS in 
Fig. 2b shows that the Ti ‘idxy is solely responsible for 
the interfacial metallic states as in the case of the well- 
known LAO/STO HS system, [35, 36] which agrees well 


with the d 2 ,j,-character-like orbital shape in the charge 
density plot, while the dyz and d^z orbitals remain un¬ 
occupied and stay at higher energies in the conduction 
band. In contrast, the orbital-resolved DOS of the sur¬ 
face CaO layer exhibits that all the O 2p orbitals cross 
the Fermi level, forming the metallic states. In short, our 
theoretical calculations for the HS model with CaO sur¬ 
face termination provide an excellent explanation for the 
recent experimental 2DEG in the nonpolar CZO/STO 
HS system. [27] 


2. ZrOi Surface Termination 

Next, we studied Ti02/Ca0 interface model with Zr02 
surface termination. We calculated the total DOS for 
unrelaxed (first column) and relaxed (second column) 
(CZO)„/STO {n = 5, 6, 7, and 8) HS models in Eig. 
3. As in the case of the HS model with CaO surface ter¬ 
mination, all the unrelaxed HS models with Zr02 surface 
termination exhibit an insulating behavior (first column 
in Eig. 3). In contrast, the relaxed HS models with n 
> 8 show a typical conducting character, with a number 
of states crossing Fermi level (second column in Fig. 3). 
As discussed later, the conducting property is induced by 
the polarization-driven charge transfer from the middle 
CaO layer in the CZO film to the interfacial Ti02 lay¬ 
ers, which lead to the partially unoccupied O 2p orbitals 
and the partially occupied Ti 3(i orbitals and thus p-type 
and n-type conducting states, that is, the co-existence 
of Two-Dimensional-Hole-Gas (2DHG) and 2DEG. For 
the (CZO)„/STO HS models with n € {5, 6, and 7}, the 
Fermi level is pinned close to the valence band maximum, 
and there exist only minor states just above the Fermi 
level, showing p-type electronic structure character. The 
partial DOS analysis for these systems shows that the mi¬ 
nor gap states arise from the partially unoccupied O 2p 
orbital and Zr 4c? orbital, see Fig. IS in the Supporting 
Information. As discussed later, this is because the CZO 
film form a “up” polarization domain toward the vacuum, 
leading to the charge transfer from the O 2p orbital of the 
STO substrate to the Zr Ad orbital of the surface Zr02 
layer. Similar electronic property also appears in the 
HS model with less CZO layers, i.e., (CZO)„/STO (n=3 
and 4), while the HS model (CZO) 2 /STO exhibits an 
insulating behavior, indicating an insulating-to-p-type- 
conductivity transition at n=3, see Fig. 2S of the Sup¬ 
porting Information. In the (CZO) 2 /STO model, as con¬ 
firmed from the structural analysis, the polarization in 
the CZO film is much weaker than that in other mod¬ 
els, which is not capable of driving the charge transfer 
between the CZO film and STO substrate and thus the 
system keeps insulating. 

The analysis for these relaxed HS structures indicates 
that the Ca and Zr cations in the CZO film move towards 
the vacuum, producing relative displacements between 
the cations and anions and the resulting “up” polariza¬ 
tion, Pup, which lead to the minor charge transfer from 
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FIG. 3. (Color online) Calculated total DO for unrelaxed 
(first column) and relaxed (second column) (CZO)rt/STO (n 
= 5, 6, 7, and 8) HS-based slab (Ti02/Ca0 interface model) 
with Zr 02 surface termination. 


the Ti02 layer to the surface Zr02 layer, see Fig. 3S 
in the Supporting Information. It is worth emphasizing 
that the HS models (CZO)„/STO (n = 3, 4, 5, 6, and 7) 
do not form 2DEG because their conductivity is caused 
by the partially unoccupied O 2p orbitals instead of Ti 
3d orbitals. Moreover, their resistivity is expected to be 
much higher than that of (CZO)„/STO (n > 8) because 
of relatively low carrier density, as shown in the DOS plot 
of Fig. 3. 

To further explore the formation mechanism of the 
2DEG in the Ti02/Ca0 interface model with Zr02 sur¬ 
face termination, we calculated the layer-resolved DOS of 
the (CZ0)8/ST0 HS system in Fig. 4, along with the 3D 
charge density projected on the bands forming the metal¬ 
lic states near the Fermi level. One can clearly see that 
the metallic states mainly come from two components: 
the partially unoccupied O 2p orbitals from the middle 
CaO layer and adjacent Zr02 layers, and the partially 
occupied Ti 3d orbitals near the interfacial region in the 
STO substrate, along with few partially occupied Zr 4d 
orbitals on the surface. Interestingly, in contrast to the 
model with GaO surface termination, the (GZO)„/STO 
HS with n > 8 model with Zr02 surface termination ex¬ 
hibits two polarization domains with opposite directions, 
as shown in the charge density plot and Fig. 3S in the 
Supporting Information. Above the middle GaO layer, 
the and Zr'^^ cations in the GZO film move to¬ 

wards the vacuum, leading to a “up” polarization, Pu-p- 
This strong polarization induces a charge transfer from 
the occupied O 2p orbitals of the middle CaO layer to 
the surface Zr 02 layer, producing n-type conductivity 
on the surface Zr 02 layer and p-type conductivity in 
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FIG. 4. (Color online) Calculated (a) layer-resolved par¬ 
tial DOS for relaxed (CZO)8/STO slab (Ti02/Ca0 interface 
model) with Zr 02 surface termination, along with the charge 
density projected on the bands forming the metallic states 
near the interface and surface region, (b) Orbital-resolved 
partial DOS for Ti atoms at the IF-I Ti02 layer and for O 
atom from the middle CaO layer. 


the middle CaO layer, see Fig. 4a and Fig. 3S of the 
Supporting Information. In contrast, below the middle 
CaO layer, the cations of the GZO film move towards the 
STO substrate, and the relative displacement between 
the cations and anions produces a “down” polarization 
Pdn that causes the charge transfer from the GZO film 
to the STO substrate. It is noted that, similar to the 
case of the HS model with CaO surface termination, the 
STO substrate exhibit a much weaker distortion within 
the first three Ti02 layers than the GZO film, and the 
deep Ti02 layers nearly show no distortion. This im¬ 
plies that a polarization discontinuity exists between the 
GZO film and STO substrate. Therefore, the transferred 
charge from the middle CaO layer to the STO substrate 
was captured by the localized Ti 3c? orbitals in the first 
few interfacial Ti02 layers of the STO substrate, forming 
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FIG. 5. (Color online) Calculated total DOS for unrelaxed 
(first column) and relaxed (second column) (CZO)n/STO (n 
= 3, 4, 5, and 6) HS-based slab (Sr0/Zr02 interface model) 
with CaO surface termination. 
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the interfacial metallic states. The partial DOS analysis 
displays that only the first three Ti02 layers, i.e., IF-I, 
IF-III, and IF-V layers, contribute to interfacial conduct¬ 
ing states, which confirms the formation of the 2DEG. 
The calculated orbital-resolved partial DOS plot shows 
that the interfacial metallic states on the IF-I layer only 
come from the dxy orbitals of Ti ions, see Fig. 4b. It is 
particularly worth mentioning that, in the CZO/STO HS 
system with Zr02 surface termination, besides the 2DEG 
in the interfacial region, there also exists a 2DHG in the 
middle layer of the GZO film, which is contributed by the 
Px/Py orbitals of O ions. This is because the “down” and 
“up” polarization in the GZO film leads to the charge 
transfer from GZO film to the STO substrate and the 
surface Zr 02 layer simultaneously, thus producing the p- 
type hole states in the middle CaO and adjacent Zr02 
layers. 
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FIG. 6. (Color online) Calculated layer-resolved partial DOS 
for relaxed (CZO) 5 /STO slab (Sr0/Zr02 interface model) 
with CaO surface termination, along with the charge den¬ 
sity projected on the bands forming the metallic states near 
the interface and surface region. 


no polarization in the GZO film. After structural relax¬ 
ation, the (CZO)„/STO models with n £ {3 and 4} ex¬ 
hibit insulating behavior, and as n increases up to 5, an 
insulator-to-metal transition occurs. This indicates that, 
for the Sr0/Zr02 interface model with the CaO surface 
termination, 5 GZO uc is the critical thickness for form¬ 
ing conductivity, which is less than that of the Ti 02 /Ca 0 
interface model with the same surface termination. 


B. Sr0/Zr02 Interface Model 

1. CaO Surface Termination 

In this section, we turn to the case of Sr0/Zr02 in¬ 
terface model by beginning with the HS model with 
CaO surface termination. To explore the possible crit¬ 
ical thickness for forming interfacial metallic states, we 
performed calculations for (CZO)„/STO (n = 3, 4, 5, 
and 6) HS model. The calculated total DOS for the un¬ 
relaxed and relaxed HS systems are shown in Fig. 5. As 
expected, before structural relaxation, all the unrelaxed 
HS systems exhibit insulating character because there is 


To investigate the origin of the metallic states in 
the Sr0/Zr02 interface model, we calculated its layer- 
resolved DOS and 3D charge density plot projected on 
the bands forming the metallic states, see Fig. 6. Similar 
to the case of Ti 02 /Ca 0 interface model, the compres¬ 
sive strain applied in the GZO film causes a downward 
movement of the Ga^"*" and Zr^+ cations to the interface, 
leading to a “down” polarization, Pdn- This Pdn leads 
to a charge transfer from the occupied O 2p orbitals in 
the surface GaO layer to the Ti 3(i orbitals at the in¬ 
terface. The transferred electrons are mainly captured 
by the localized Ti 3d orbitals of the first three Ti02 
layer in the STO substrate, producing n-type interface 
conducting states. 
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FIG. 7. (Color online) Calculated total DOS at for unrelaxed 
(first column) and relaxed (second column) (CZO)rt/STO (n 
= 8.5, 9.5, 10.5, and 11.5) HS-based slab (Sr0/Zr02 interface 
model) with Zr 02 surface termination. 


2. ZrOi Surface Termination 


Next, we studied the Sr0/Zr02 interface model with 
Zr02 surface termination. The calculated total DOS for 
the unrelaxed and relaxed (CZO)„/STO (n = 8.5, 9.5. 
10.5, and 11.5) HS model is shown in Fig. 7. Before 
the structural relaxation, all the models show an insulat¬ 
ing property. After structural relaxation, the CZO/STO 
HS models with n = 8.5 and 9.5 show a p-type conduc¬ 
tivity, in which the Fermi level is pinned slightly below 
the valence band maximum with minor states crossing 
the Fermi level. Similar electronic structure character 
is also found in the (CZO)„/STO HS models with n 
G {4.5, 5.5, 6.5, and 7.5}. The partial DOS analysis 
informs that the minor states above the Fermi level in 
these HS models come from the unoccupied O 2p orbital 
in the STO substrate. This is because, in these relaxed 
(CZO)„/STO HS models (n = 4.5 to 9.5), the“up” po¬ 
larization, Pup-, caused by the relative displacements be¬ 
tween the cations and anions, leads to a charge transfer 
from the STO substrate to the surface Zr02 layer (see 
Fig. 4S in the Supporting Information), similar to the 
case of the Ti 02 /Ca 0 interface model with Zr 02 surface 
termination (Fig. 3S in the Supporting Information). 

At n > 10.5, the (CZO)„/STO HS system exhibits 
metallic states. The calculated layer-resolved DOS of 
the (CZO)io. 5 /STO HS model and 3D charge density 
projected on the bands forming the metallic states are 
shown in Fig. 8. It shows that two polarization domains 
with opposite directions and Pdn in the CZO film, 
separated by the middle CZO layer, are formed in the 


ZrO^ surface termination 



FIG. 8. (Golor online) Calculated layer-resolved partial DOS 
for relaxed (CZO)io. 5 /STO slab (Sr0/Zr02 interface model) 
with Zr 02 surface termination, along with the charge density 
projected on the bands forming the metallic states near the 
interface and surface region. 


CZO film. In the Pup domain, the Ca^"*" and Zr^+ cations 
move towards the vacuum, and the polarization induces 
a charge transfer from the middle CaO layer to the sur¬ 
face Zr 02 layer, creating n-type conductivity on the sur¬ 
face Zr 02 layer and p-type conductivity in the middle 
CaO layer, see Fig. 8. In the Pdn domain, in contrast, 
the cations of the CZO film move towards the STO sub¬ 
strate, and the Pdn polarization leads to a charge transfer 
from the middle CaO layer to the STO substrate, produc¬ 
ing u-type interfacial conductivity. This indicates that, 
in the Sr0/Zr02 interface model with Zr02 surface ter¬ 
mination, the two polarization domains lead to the co¬ 
existence of the 2DEG on the interface and the 2DHG 









































































0.8 




no. of CaZrOj unit cells 


FIG. 9. (Color online) Calculated average relative displace¬ 
ment (Sz) of Ca and Zr cations with respect to O ions in the 
corresponding CaO and Zr02 planes in the CZO film for the 
Ti02/CaO interface model with (a) CaO and (b) Zr02 surface 
termination. Vertical dashed green lines indicate the critical 
thickness of the CZO films for forming 2DEG at the interface 
in each model. For the (CZO)n/STO HS model with Zr02 
surface termination, at n> 8, two polarization domains, Pup 
and Pdn, are formed, and the Ca/Zr cations move in opposite 
directions in the two domains. 


in the CZO film. The critical thickness of the CZO film 
for forming such an electronic state is about 10.5 uc. In 
short, our results indicate that over all, the polarization 
and electronic structure characters of the Sr0/Zr02 in¬ 
terface model are similar to those of the Ti 02 /Ca 0 inter¬ 
face model for the each type of the surface termination, 
though the required critical thickness to form the inter¬ 
facial conducting states in the STO substrate is different. 
Hence, it is expected that the 2DEG can also be observed 
in the Sr0/Zr02 interface model, as in the case of the 
Ti02/Ca0 interface model. [27] 


C. Structural Relaxation 

Since the Ti02/Ca0 interface model shares similar 
polarization and electronic structure characters with 
those of the Sr0/Zr02 interface model, in the follow¬ 
ing discussions, we only take the Ti 02 /Ca 0 interface 
model as an example to explore the polarization prop¬ 
erty and electron carrier density for CZO/STO HS sys¬ 
tem. As discussed above, the structural optimization 


plays a crucial role in producing the metallic states in 
all the CZO/STO HS models. Specifically, the strain- 
induced relative displacement between the cations and 
anions, i.e., polarization, leads to the formation of 2DEG 
(2DHG) in the CZO/STO HS system. To quantitatively 
characterize the relationship between the polarization 
and the polarization-induced 2DEG (2DHG) in the two 
CZO/STO HS models, we plotted the average relative 
displacement {5z) of the Ca/Zr cations in the CZO film 
along the c-axis as a function of the number of the CZO 
uc (n) in Fig. 9 in each case. The Sz of Qj&jJx cations 
was calculated with respect to the oxygen ions in the cor¬ 
responding CaO and Zr02 plane. The positive and nega¬ 
tive values indicate that Ca/Zr cations move towards and 
away from the interface, respectively. Our results clearly 
show that, in the HS model with CaO surface termina¬ 
tion, Szca > 0 for the whole range of n, indicating that 
Ca cations move towards the interface, while Zr cations 
slightly move away from the interface at n < 6.5, see 
Fig. 9a. This implies that the total polarization in the 
CZO film is almost negligible below the critical thickness. 
When n increases up to 6.5, be., the critical thickness of 
the CZO film for forming the 2DEG, Sz is about 0.65 A, 
and as n further increases, Sz slightly increases. For the 
(CZO)„/STO HS model with Zr02 surface termination, 
as mentioned above, at n < 7, the Ca/Zr cations move 
away from the interface, be., Sz < 0. As shown in Fig. 9b, 
when n increases from 2 to 7, j i5z j gradually increases 
up to 0.5 A, and as n further increases, j Sz j nearly keeps 
constant. At n > 8, two polarization domains, Pup and 
Pdn, are formed, as shown in Fig. 4a. That is to say, 
the Ca/Zr cations move from the middle CZO layers to¬ 
wards the vacuum and STO substrate, respectively, in 
the Pup and Pdn polarization domain, and | (5z | in the 
two polarization domains is comparable. 

To quantitatively characterize the polarization 
strength, we calculated the polarization P of the 
CZO film for the two HS models using the following 
formula: [37, 38] 


N 

p = ( 1 ) 

^ i=l 

where H is the volume of the CZO film, N is the 
number of cation in the CZO film, Z* is Born effective 
charge for each cation (Ca/Zr), and Szi is the relative dis¬ 
placement between the ith cation (Ca/Zr) and the anion 
(O) in the corresponding CaO/Zr02 plane, be., Szi=zo- 
zca/Zr- Our Calculated Born effective charges Z* are 
2.27 and 4.93 for Ca and Zr, respectively, and —3.88 and 
— 1.66 for O ions in the CaO and Zr02 planes, respec¬ 
tively, for the tetragonal bulk CZO. For the CZO/STO 
HS model with CaO surface termination, at n=6.5, be., 
critical thickness of forming 2DEG, the estimated polar¬ 
ization is about 66.5 /rC/cm^. For the HS model with 
Zr 02 surface termination, at ?r=8, two polarization do¬ 
mains are formed in the CZO film. The polarization of 
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each domain was calculated as about 57.4 ^C/cm^ for 
Pdn and 46.4 fiG/cm^ for P^p, respectively. These values 
are comparable with the the closely related SrZrOa/STO 
superlattice system, in which the polarization was esti¬ 
mated to be around 41-43 ^C/cm^.[39, 40] Nevertheless, 
it is noted that the polarization value calculated from 
first-principles calculation is much larger than the exper¬ 
imental value of 3.5 /xC/cm^.[27] This is probably be¬ 
cause that, in the epitaxial materials growth process, as 
the thickness of the CZO film increases, the CZO film un¬ 
dergoes less compressive strain from the STO substrate 
and tends to recover its bulk equilibrium condition, and 
thus the relative displacement 6z of the cations in the 
CZO film decreases significantly. In contrast, in the com¬ 
putational simulation, the lattice constant of the whole 
HS model in the a6-plane was fixed due to the period¬ 
icity requirement in the plane-wave-based DFT calcula¬ 
tions, which means that the whole CZO film undergoes 
the same compressive strain due to the lattice mismatch, 
thus leading to the larger average 5z than the experimen¬ 
tal value. Hence, this discrepancy is expected to become 
larger as the CZO film grows thicker. 

In short, our calculations show that, when the CZO 
film exceeds the critical thickness, i.e., 6.5 uc. for the 
model with CaO surface termination and 8 uc. for the 
model with Zr02 surface termination, both the two HS 
systems can lead to polarization, Pdn, and the result¬ 
ing 2DEG on the interface. The difference is that only 
one polarization domain, Pdn, is formed in the model 
with CaO surface termination, while two polarization 
domains, i.e., Pup and Pdn, are formed in the model 
with Zr02 surface termination. This is probably be¬ 
cause the CaO-terminated surface and the interface share 
the same role in driving the Pdn towards the STO sub¬ 
strate, while the Zr02-terminated surface and the inter¬ 
face have a competing role in driving the polarization 
due to the electrostatic potential discontinuity. In ad¬ 
dition, the CZO film with CaO surface termination was 
predicted to be energetically more stable than that with 
Zr02 surface termination, [41] implying that the experi¬ 
mental CZO/STO HS sample is likely to be terminated 
with CaO surface. In fact, as discussed above, our pre¬ 
dicted 2DEG in the HS-based slab model with CaO sur¬ 
face termination is more consistent with the Chen et a/.’s 
experimental observation. [27] 


D. Interfacial Electron Carrier Density 

To examine the influence of the CZO film thickness 
on the interfacial charge carrier density, we computed 
the occupation number of the Ti 3d orbitals from the 
interfacial (lE-I) Ti02 layer by integrating the partial 
DOS of the Ti 3d occupied states near the Eermi level 
as well as the interfacial charge carrier density for the 
Ti02/Ca0 interface model with CaO and Zr02 surface 
terminations in Fig. 10. This is because that the IF-I 
Ti 02 layer plays a major role in the interfacial conductiv- 
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FIG. 10. (Color online) Calculated orbital occupation and 
electron carrier density of Ti 3d orbitals from the IF-I Ti02 
layer in the STO substrate with respect to the number of 
CZO unit cells (n) for the Ti02/CaO interface model with 
(a) CaO and (b) Zr02 surface termination. Vertical dashed 
green lines indicate the critical thickness of the CZO films for 
forming 2DEG at the interface in each model. 


ity in the CZO/STO HS system. For the (CZO)„/STO 
HS model {n = 3.5, 4.5, 5.5, 6.5, 7.5, 8.5, and 9.5) with 
CaO surface termination, the calculated orbital occupa¬ 
tion of Ti 3d orbital and interfacial charge carrier density 
as a function of the number of CZO uc {n) was plotted 
in Fig. 10a. Similarly, for the (CZO)„/STO HS model {n 
= 4, 5, 6, 7, 8, 9, and 10) with Zr02 surface termination, 
the corresponding orbital occupation number and inter¬ 
facial charge carrier density are shown in Fig. 10b. As 
discussed above, below the critical thickness of the CZO 
film for forming the 2DEG at the interface, there is no 
electrons captured by the Ti 3d orbitals in the IF-I Ti02 
layer, and thus the systems exhibit zero orbital occupa¬ 
tion and interfacial electron carrier density. Above the 
critical thickness of the CZO film, the interface becomes 
conductive, and as the CZO unit cell n increases, their 
interfacial electron carrier density increases., which is in a 
good agreement with the recent experimental finding. [27] 


IV. CONCLUSION 


In summary, we studied the formation of 2DEG 
in the CZO/STO HS-based slab systems using first- 
principles electronic structure calculations. Four types of 
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CZO/STO HS-based slab systems, including Ti02/Ca0 
and Sr0/Zr02 interface models with CaO and Zr02 sur¬ 
face terminations in each model, were modelled. It was 
found that, as the thickness of CZO film increases, the 
compressive strain exerted on the CZO film, caused by 
the lattice-mismatch between the CZO film and STO 
substrate, leads to a strong polarization in the CZO layer 
and the resulting interfacial 2DEG at the interface. Nev¬ 
ertheless, the polarization direction and critical thickness 
of the CZO film for forming the 2DEG strongly depend 
on the surface termination of the CZO film. For the 
HS model with CaO surface termination, a critical thick¬ 
ness of CZO film of approximately 6.5 and 5 unit cells 
is necessary to produce a strong polarization towards the 
STO substrate in the Ti02/Ca0 and Sr0/Zr02 interface 
models, respectively. The polarization induces the charge 
transfer from the surface CaO layer to the hrst few Ti02 
layers in the STO substrate, thus forming the 2DEG at 
the interface. For the HS model with Zr02 surface termi¬ 
nation, in contrast, as the CZO film increases up to 8 unit 
cells for Ti02/Ca0 interface model and 10.5 unit cells for 
Sr0/Zr02 interface model, the lattice-mismatch-induced 


strain produces two polarization domains with opposite 
directions in the CZO film, which points towards the in¬ 
terface and surface from the middle CZO layer, respec¬ 
tively. The two polarization domains leads to the co¬ 
existence of the 2DEG at the interface and the 2DHG 
at the middle CZO film, along with the n-type few sur¬ 
face metallic states. Our results reveal the mechanism of 
producing the high-mobility 2DEG (2DHG) in the non¬ 
polar CZO/STO HS system based on the stain-induced 
polarization discontinuity, which provides a new avenue 
to produce 2DEG (2DHG) in the perovskite-based HS 
systems. 
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